ACCEPTED TO Astrophysical Journal Letters 

Preprint typeset using IAT^X style cmulatcapj v. 4/12/04 



THE GALEX-VVDS* MEASUREMENT OF THE EVOLUTION OF THE 1500A LUMINOSITY FUNCTION 



S. Arnouts \ D. Schiminovich 2,3 , O. Ilbert 1 , L. Tresse l , B. Milliard \ M. Treyer 

10 



S. Bardelli 



T. 



^1" 

o 
o 

(N 
> 

o 
in 



o 

Or 

o 

CO 

c3 



Budavari , T. K. Wyder 3 , E. Zucca , O. Le Fevre , C. Martin , G. Vettolani , C. Adami , M. Arnaboldi 
15 , T. Barlow 3 , L. Bianchi 4 , M. Bolzonella ls , D. Bottini n , Y.-I. Byun 5 , A. Cappi 17 , S. Charlot 13,16 , T. 
Contini 12 , J. Donas 1 , K. Forster 3 , S. Foucaud n , P. Franzetti 11 , P. G. Friedman 3 , B. Garilli 11 , I. 



Gavignaud 

D. Maccagni 11 , B. F. Madore ', R. Malina i , B. Marano ls , C Marinoni \ H.J. McCracken iB , A. Mazure \ B 



L. Guzzo , T. M. Heckman , C. Hoopes , A. Iovino 

18 



P. Jelinsky °, V. Le Brun \ Y.-W. Lee 

16 



Meneux 1 , 
pozzetti 1t , 



R. Merighi , P. Morrissey , S. Neff s , S. Paltani , R. Pello , J. P. Picat , A. Pollo 
M. Radovich I6 , R. M. Rich 9 , R. Scaramella io , M. Scodeggio n , M. Seibert 3 , O. Siegmund b , T. 
Small 3 , A. S. Szalay 



B. Welsh 8 , C. K. Xu 19 , G. Zamorani 17 , A. Zanichelli iu 

accepted to Astrophysical Journal Letters 

ABSTRACT 

We present the first measurement of the galaxy luminosity function at 1500 A between 0.2 < z < 1.2 
based on GALEX-VVDS observations (~1000 spectroscopic redshifts for galaxies with NUV < 24.5) 
and at higher z using existing datasets. Our main results are summarized as follows : (i) luminosity 
evolution is observed with AM* ~ —2.0 mag between z — and z = 1 and AM* ~ —1.0 mag between 
z = 1 and z = 3. This confirms that the star formation activity was significantly higher in the past, 
(ii) the LF slopes vary between —1.2 > a > —1.65, with a marginally significant hint of increase 
at higher z. (iii) we split the sample in three restframe (B — I) intervals providing an approximate 
spectral type classification: Sb-Sd, Sd-Irr and unobscured starbursts. We find that the bluest class 
evolves less strongly in luminosity than the two other classes. On the other hand their number density 
increases sharply with z (~15% in the local universe to ^55% at z ~ 1) while that of the reddest 
classes decreases. 

Subject headings: galaxies: evolution — galaxies: luminosity function — ultraviolet: galaxies - 
cosmology: observations 



*THE VVDS OBSERVATIONS HAVE BEEN OBTAINED WITH 
THE EUROPEAN SOUTHERN OBSERVATORY VERY LARGE 
TELESCOPE, PARANAL CHILE 
Electronic address: stephane.arnouts@oamp.fr 

Laboratoire d'Astrophysique de Marseille, BP 8, Traverse du 
Siphon, 13376 Marseille Cedex 12, France 

2 Department of Astronomy, Columbia University, MC2457, 550 
W. 120 St. NY,NY 10027 

3 California Institute of Technology, MC 405-47, 1200 East Cal- 
ifornia Boulevard, Pasadena, CA 91125 

4 Department of Physics and Astronomy, The Johns Hopkins 
University, Homewood Campus, Baltimore, MD 21218 

5 Center for Space Astrophysics, Yonsei University, Seoul 120- 
749, Korea 

6 Space Sciences Laboratory, University of California at Berke- 
ley, 601 Campbell Hall, Berkeley, CA 94720 

7 Observatories of the Carnegie Institution of Washington, 813 
Santa Barbara St., Pasadena, CA 91101 

8 Laboratory for Astronomy and Solar Physics, NASA Goddard 
Space Flight Center, Greenbelt, MD 20771 

9 Department of Physics and Astronomy, University of Califor- 
nia, Los Angeles, CA 90095 

10 INAF-IRA, Via Gobetti,101, 1-40129, Bologna, Italy 

11 INAF-IASF, via Bassini 15, 1-20133, Milano, Italy 

12 Laboratoire d'Astrophysique de l'Observatoire Midi-Pyrenees, 
14, avenue E. Belin, F-31400 Toulouse, France 

13 Max Planck Institut fur Astrophysik, 85741 Garching 
many 

14 INAF-Osservatorio Astronomico di Brera, via Brera 28,. 
Ian, Italy 

15 INAF-Osservatorio Astronomico di Capodimontc, 
Moiariello 16, 1-80131 Napoli, Italy 

16 Institut d'Astrophysique de Paris, UMR 7095, 98 bis 
Arago, F- 75014 Paris, France 

17 INAF-Osservatorio Astronomico di Bologna, via R,anzani,l, 
1-40127 Bologna, Italy 

18 Universita di Bologna, Departimento di Astronomia, via Ran- 
zani 1, 1-40127 Bologna, Italy 

19 NASA/IPAC, California Institute of Technology, Mail Code 
100-22, 770 S. Wilson Ave., Pasadena, CA 91125 



Ger- 



Mi- 



bvd 



1. INTRODUCTION 

Evidence collected over the past decade suggests that 
the cosmic, volume-averaged star formation rate (SFR) 
density has increased from the present to z ~ 1 (Lilly 
et al. 1996, Wilson et al. 2002) with a possible flatten- 
ing at higher redshift (Steidel et al. 1999, Giavalisco et al 
2004). The details of the evolution between < z < 1 are 
still being debated and appear to depend on the selected 
wavelength. Restframe far-ultraviolet (FUV: 1500A) ob- 
servations can provide a sensitive measurement of the 
ongoing SFR but up to now have only been obtained in 
the local universe (2000A, Sullivan et al. 2000) and at 
z > 2.5 (1700A, Steidel et al. 1999). The GALEX mis- 
sion now allows us to refine the local FUV measurements 
(Wyder et al. 2004, Treyer et al. 2004, Budavari et al. 
2004) and to fill the redshift gap where most of the SFR 
evolution has taken place. 

In this paper we use a unique spectroscopic sample of 
^1000 NUV-selected galaxies to derive the first measure- 
ment of the 1500A LF between 0.2 < z < 1.2. We also 
use HDF data to extend our analysis to z = 3. These 
results allow us to address the evolution of the lumi- 
nosity density and SFR, presented in a companion pa- 
per (Schiminovich et al. 2004). Throughout the paper, 
we assume a flat ACDM cosmology with J7m = 0.3 and 
H = 70 km s- 1 Mpc" 1 . 

2. DATA DESCRIPTION 

We use the data collected in the 02hr field 
(02h26m00s,-04°30'00") by GALEX and by the VIR- 
MOS VLT Deep Survey (VVDS). The 02hr field was 
observed as part of the GALEX Deep Imaging Survey 
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(DIS), with a total integration time of t exp = 52763sec 
in two channels (FUV -1530A and NUV ~2310A). We 
restrict the analysis to the central 1° diameter region 
of the GALEX field of view which provides higher uni- 
formity (Morrissey et al. 2004). The typical PSF has 
a FWHM ~ 5 arcsec. Because of the high source 
density in these images, the GALEX data pipeline (im- 
plementing a modified version of SExtractor, Bertin & 
Arnouts 1996) was run with deblending parameters con- 
figured to improve extraction of unresolved point sources 
at the faint limit. The number counts are ~80% com- 
plete at 24.5 (AB system) in the FUV and NUV bands 
(Xu et al. 2004). Magnitudes were corrected for Galac- 
tic extinction using the Schlegel et al. (1998) redden- 
ing map and Cardelli et al. (1989) extinction law with 
A FUV /E(B -V) = 8.29, A NU y/E(B -V) = 8.61. 

The VVDS is a spectroscopic survey with multicolor 
photometry (Le Fevre et al 2003). Deep B,V,R,I 
imaging over 1.2 deg 2 was carried out using the 
Canada-France-Hawaii telescope, with typical depths of 
27,26,26.2,25.3 respectively and a PSF better than 1 arc- 
sec (McCracken et al. 2003). The ongoing deep spectro- 
scopic survey covers ^1950 arcmin 2 with a preliminary 
sample of ^7200 spectroscopic redshifts of galaxies with 
17.5 < I A b < 24. 

In the UV sample, all sources have at least one optical 
counterpart (hereafter OC) within a radius of 4". With 
the VVDS spectroscopic sample, we match 1157 galax- 
ies with NUV < 24.5 and find that 48% are isolated 
(one OC in 4"), 36% have two OCs and 16% have more 
than two OCs. These fractions are almost constant with 
magnitude between 21.5 < NUV < 24.5 and with red- 
shift. As a preliminary matching procedure, we associate 
the closest optical counterpart to each UV source. This 
also selects the brightest B band counterpart in 95% of 
the double counterpart cases and 85% of the more than 
two counterpart cases. In the case of NUV sources with 
several OCs, the UV flux may result from the UV con- 
tribution of several sources. Although we do not correct 
for the UV flux attributed to the chosen OC, we use sim- 
ulations to estimate the magnitude of the bias. 

To test the impact of blends, we divide the NUV flux 
among the potential OCs determined using a match re- 
liability statistic based on the likelihood estimator of 
Sutherland and Saunders (1992). We find that the aver- 
age UV flux overestimate is ~ —0.25 mag in the case of 
two OCs and up to ~ —0.5 mag for the more complex 
cases. By comparing the luminosity functions derived 
for this blend-corrected NUV selected sample and the 
original sample, we find that the variations of is less 
than —0.25, within the errorbars quoted below. This 
method may be more accurate but still requires several 
strong assumptions (optical and UV sources are intrinsi- 
cally unresolved, optical- UV color is fixed). As we have 
shown that the results are not significantly affected by 
the adopted method, we use the simplest approach. 

Figure ^ shows the (NUV — I) vs NUV color magni- 
tude diagram. Isolated objects are displayed with filled 
circles and multiple counterpart cases with crosses. The 
spectroscopic sources (red symbols) randomly sample the 
full range of colors observed in the total sample and is 
therefore unbiased. For / band saturated galaxies (~4% 
with Iab < 17.5, green symbols), we estimate that ~95% 
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Fig. 1. — Color-magnitude diagram {(NUV-I)ab vs NUVab)- 



are at z < 0.2 using photometric redshifts in GALEX- 
SDSS fields (Budavari et al. 2004). No redshifts were 
measured for galaxies with Iab > 24, selecting against 
very blue objects, but the small fraction (~4%) is un- 
likely to affect our results. 

Although our preliminary catalog contains galaxies 
with spectroscopic redshifts between < z < 1.5 we 
restrict the sample to the 1039 galaxies in the redshift 
range between 0.2 < z < 1.2. The low redshift cut avoids 
local galaxies saturated in the / band. The high redshift 
limit ensures that the the 912A break falls below the blue 
edge of the NUV passband. 

3. THE FUV (1500A) LUMINOSITY FUNCTIONS 

We derived restframe FUV absolute magnitudes using 
k-corrected NUV fluxes (with k-correction based on SED 
fits). This choice yields minimal k-corrections because 
the average wavelength of the NUV passband samples 
the rest-FUV interval 1050A< X rest < 1925A over our 
adopted redshift range. To measure the luminosity func- 
tion (LF) we use the VVDS LF tool (ALF) by Ilbert et al. 
2004a, which includes the l/V max , C+, SWML and STY 
estimators. A weight is assigned to each spectroscopic 
galaxy to take into account the spectroscopic strategy of 
the VVDS observations and the completeness correction 
of the NUV number counts. We refer to a forthcoming 
paper for a complete description of the absolute mag- 
nitude and weight measurements. Figure shows the 
GALEX- VVDS LFs in four redshift bins (4 upper pan- 
els) for the four estimators (V max : circles, SWML: trian- 
gles, C+: squares, STY: solid line). For the STY fits, we 
show the extrema (hashed area) based on the individual 
la uncertainty in a and M+ (the likelihood probability 
contours are shown as insets for 2AlnC = 1: dashed 
lines and for 2AlnC — 2.3: solid lines). The vertical 
lines show the limits beyond which the global sample be- 
comes biased against certain spectral types (see Ilbert et 
al. 2004b). These values are listed in Table \\\&s Muas- 
Points beyond these limits are not included in the fitting 
procedures. In general, NUV selection serves to reduce 
this bias significantly. Finally we measure the FUV-LF 




Redshift evolution of the 1500A luminosity function 
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Fig. 2.— 1500A LF between 0.2 < z < 3.5. 



at high-z using photometric redshifts in the HDF North 
and South (Arnouts et al. 2002; bottom panels). For the 
purpose of restframe FUV selection, we define one sam- 
ple between 1.75 < z < 2.25 using F450 passband with 
F450as < 27 and one sample between 2.40 < z < 3.40 
using F606 passband with F606ab < 27. At z ~ 3, we 
also compare our results with the LF at 1700A (Steidel 
et al. 1999, red dashed line). As a reference, we show 
the local 1500 A LF derived from GALEX data (Wyder 
at al. 2004; dotted lines). All of the STY parameters 
are listed in Table Q] and Figure 03 shows the slope (a, 
top panel) and the M* (bottom panel) parameters ver- 
sus redshift (GALEX- VVDS sample: filled circles, HDF 
sample: stars, LBG sample: open squares, local sample: 
filled squares). 

The faint end slopes vary between —1.65 < a < —1.2 
for 0.2 < z < 3 with a marginal steepening with z 
(within the la errorbars). The M* vs z plot reveals 
strong redshift evolution. A significant brightening of 
order AM* ~ —2 occurs between < z < 1.2. The 
higher z samples show that the trend continues to z ~ 3 
at a lower rate (AM* ~ —1). This increase is highly 
significant with respect to our errorbars and the source 
blending issue discussed in Section|2| The brightening of 
M* and the steepening of the slope observed at 1500 A 
is qualitatively consistent with the evolution detected at 
longer wavelengthes (2800A, U and B bands) by Wolf et 
al. (2003) and Ilbert et al. (2004a). 

4. DEPENDENCE OF THE FUV LF ON COLOR 

The luminosity function of galaxies has been shown to 
vary as a function of rest-frame color or spectral type 
(Blanton et al. 2001, Wolf et al. 2003). In this section 
we explore the FUV LF and its evolution as a function 
of color. 

As shown by Salim et al. (2004), restframe (NUV-R) 



GALEX-AIS A Red type « HDFs 

CALEX-VVDS O Sd-Irr type □ LBCs 
□ SB type 



i — r 



Fig. 3. — Evolution of the LF's STY parameters vs redshift. 



TABLE 1 

1500A LF PARAMETERS FOR (firj = 0.3, Qa 
KM S _1 MPC -1 



: 0.7) and H = 70 
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0.055" 
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-1.21 ±0.07 


-18.05 ± 0.11 


4.07 ±0.56 


0.2 


- 0.4 


319 


-15.36 


-1.19 ±0.15 


-18.38 ± 0.25 


6.15 ± 1.76 


0.4 


- 0.6 


258 


-16.85 


-1.55 ±0.21 


-19.49 ±0.37 


1.69 ±0.88 


0.6 


- 0.8 


274 


-17.91 


-1.60 ±0.26 


-19.84 ±0.40 


1.67 ±0.95 


0.8 
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188 
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196 




-1.50 ±0.20 
-1.50 ±0.20 


-19.531^ 
-19.63+027 
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-1.50 ±0.20 


-19.89_°i? 
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1.75-2.25-f 


139 


-17.54 


-1.49 ±0.24 


-20.33 ± 0.50 


2.65 ±2.00 


2.40-3.40-f 


173 


-18.17 


-1.47 ±0.21 


-21.08 ± 0.45 


1.62 ±0.90 


2.50-3.50 s 


564 




-1.60 ±0.13 


-21.07 ±0.15 


1.40 



10 -3 Mpc" 3 

Local sample with z < 0.1 and FUV < 20 
(B - I)> 0.85 sample with NUV < 24.5 
0.56 <(B-I)< 0.85 sample with NUV < 24.5 
(B - I) < 0.56 sample with NUV < 24.5 
HDF sample with F450 < 27 and F606 < 27 
LBG sample with R < 25 



is tightly correlated with the ratio of current to past av- 
eraged SFR, implying that the star formation history of 
a galaxy can already be constrained using a single color. 
In Fig. 01 we plot (NUV — R) vs (B — I) restframe as well 
as the location of theoretical SEDs from Elliptical to Sd 
(Poggianti 1997: filled circles), observed Irregular (Cole- 
man et al. 1980: filled circle) and Starbursts (Kinney et 
al. 1996: filled triangles). 

As (B - /) correlates with (NUV - R) (Fig gjl and 
is not affected by the UV flux uncertainties mentioned 
above, we use (B — I) as a proxy for spectral types. We 
note however the degeneracy between dust and age of the 
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Fig. 4. — Restframe color-color plot for the color classifica- 
tion (upper panel). Fraction of galaxy classes vs redshift (bottom 
panel). 



stellar population: dusty starbursts (SB3 to SB6, filled 
triangles) cannot be distinguished from spiral galaxies 
(filled circles) using (B — I). Eliminating this degeneracy 
will require additional color information sensitive to the 
old star population (e.g. K band flux) not available in 
the present dataset. 

We split our sample into three classes: (B — I) < 0.56 
(corresponding to unobscured starbursts), 0.56 < (B — 
I) < 0.85 (corresponding to Irr to Sd with a contamina- 
tion of obscured starbursts) and (B — I) > 0.85 (corre- 
sponding to Sd to Sb with a contamination of obscured 
starbursts, and a negligible 7% contribution of types ear- 
lier than Sb) . The shortage of red systems in our sample 
is consistent with the morphological analysis of de Mello 
et al. (2004) based on a small NUV restframe sample (34 
objects) in the CDF-South observed with HST (ACS). 
The authors only found 2 early type galaxies (8%), the 
rest being late types (32%) and starbursts (60%). 

In the bottom panel of Fig. we show the rela- 
tive fraction for each color class in three redshift bins: 
0.2 < z < 0.5, 0.5 < z < 0.8, 0.8 < z < 1.2. We extend 
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our analysis to lower z by applying the same criteria to 
the GALEX-SDSS fields for galaxies with z phot < 0.2 
and FUV < 22 (filled squares). We find that the con- 
tribution of unobscured starbursts rises from ~ 12% at 
z ~ 0.1 to 55% at z ~ 1, while the fractions of the two 
other classes decrease. (The reddest class actually seems 
to increase again at z ~ 1, but this may reflect a clas- 
sification bias since the / band probes the emitted light 
below the 4000A break and does not trace reliably the 
old stellar population). 

We measure the LFs for the three classes. Due to 
the small number of galaxies we fixed the slopes to 
a = —1.4 ± 0.2 for our reddest class and a = —1.5 ± 0.2 
for the two others, in agreement with the values obtained 
between 0.2 < z < 0.5. These values are consistent with 
the slopes derived by Wolf et al. (2003) at 2800A for their 
types 2 to 4. The best fit STY parameters are listed in 
Table [l]and in Fig [2 

We find that the bluest class (unobscured starbursts) 
evolves less strongly in luminosity than the two other 
classes (AM* < —0.5 mag between 0.35 < z < 1.0 as op- 
posed to AM* > — 1 mag) . On the other hand the num- 
ber densities (p = f™™* <S>(M)dM with M cut = -18.5) 
evolve similarly to the non-volume corrected fractions 
shown in Fig 0] namely the density of unobscured star- 
bursts increases sharply with z while that of the reddest 
classes decreases. 

The interpretation of our results in terms of star for- 
mation rate evolution is presented in a companion paper 
by Schiminovich et al. (2004). 
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